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SUMMARY 


The  Coastal  Engineering  Research  Center  (CERC),  under  sponsorship  from 
the  Naval  Civil  Engineering  Laboratory  (NCEL) ,  collected  wave  data  for  25  days 
from  a  closely  spaced  array  of  five  wave  gages  moored  in  17  m  of  water  off¬ 
shore  of  the  North  Carolina  coast.  The  objective  of  the  collection  was  to 
determine  the  spatial  variability  of  the  wave  field  under  a  variety  of  inci¬ 
dent  wave  conditions  to  aid  the  NCEL  in  planning  ship  mooring  tests.  The  spa¬ 
tial  variability  of  the  wave  field  was  examined  in  terms  of  the  coherence 
function  measured  between  pairs  of  wave  gages.  The  primary  factor  affecting 
the  coherence  function  is  the  spatial  separation  between  the  gages.  Coherence 
was  found  to  rapidly  decrease  as  spatial  separation  increased,  particularly  in 
the  along-crest  direction  (right-angle  to  the  direction  of  wave  travel).  The 
effect  of  spatial  separation  was  quantified  as  nondimens ional  coherence  con¬ 
tours  which  provide  a  first  rough  estimate  of  the  wave  coherence  which  would 
be  expected  between  two  gages .  Because  of  the  approximate  nature  of  the  esti¬ 
mation  technique,  collection  of  onsite  coherence  data  is  recommended  for  fu¬ 
ture  field  exercises  if  such  data  are  needed  for  subsequent  analyses.  As  a 
result  of  this  research,  some  operational  recommendations  are  provided  to  aid 
in  conducting  future  ship  mooring  experiments.  These  recommendations  are  as 
follows: 

-  Align  the  two  wave  gages  parallel  with  the  central  incident  wave 
direction. 


-  Maintain  equal  scope  on  all  wave  gage  moorings. 

-  Relocate  wave  buoys  when  possible  to  maintain  the  desired  alignment 
with  the  wave  field. 

-  Place  additional  vacant  buoy  moorings  in  an  array  to  facilitate  buoy 
relocation. 

-  Moor  a  backup  wave  buoy  at  the  experimental  site. 


□  □ 
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CERC,  Mr.  R.  P.  Savage,  former  Chief,  Research  Division,  Dr.  J.  R.  Houston, 
Chief,  Research  Division,  Dr.  C.  L.  Vincent,  former  Chief,  Coastal  Oceanogra¬ 
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SPATIAL  VARIABILITY  IN  THE  NEARSHORE  WAVEFIELD 


PART  I :  INTRODUCTION 
Background 

1.  The  Naval  Civil  Engineering  Laboratory  (NCEL)  is  investigating  ship 
motions  and  mooring  forces  for  vessels  anchored  in  coastal  waters.  A  mooring 
force  model  has  been  developed  by  NCEL,  but  field  data  are  required  in  order 
that  NCEL  can  verify  the  model's  accuracy  and  determine  its  limitations.  In 
the  NCEL  model,  the  wave  field  is  to  be  applied  at  the  center  of  gravity  of 
the  ship.  Since  this  is  an  artificial  constraint,  i.e.,  there  are  no  waves  at 
this  point,  it  is  desirable  to  estimate  the  differences  between  the  waves 
causing  ship  motion  and  those  a  few  hundred  feet  away  where  wave  data  can  be 
collected  without  excessive  ship-induced  interference.  Naturally  occurring 
waves  are  not  long-crested  and  the  waves  change  in  form  as  they  propagate. 
Since  the  form  of  each  wave  striking  the  ship  is  important,  field  studies  for 
verifying  the  mooring  force  model  could  be  planned  more  effectively  if  the 
following  information  were  available:  (a)  the  best  position  to  locate  a  gage 
to  record  wave  data  suitable  for  evaluating  mooring  response,  and  (b)  the  va¬ 
riation  in  wave  form  that  can  be  expected  between  the  recorded  waves  and  those 
which  actually  cause  the  ship  motion. 

2.  To  help  meet  this  need  for  information,  the  Coastal  Engineering  Re¬ 
search  Center  (CERC)  collected  data  for  25  days  at  the  Field  Research  Facility 
(FRF)  at  Duck,  North  Carolina,  under  NCEL  sponsorship.  Five  Waverider  wave 
measuring  buoys  were  moored  in  proximity  to  each  other,  and  wave  data  were 
collected  simultaneously  from  all  gages. 

Objective  of  the  Experiment 

3.  The  primary  objective  was  to  obtain  time  series  records  of  the  sea 
surface  displacements  for  periods  of  wind-generated  seas,  swell,  and  for  the 
case  when  more  than  one  distinct  wavetrain  was  present.  Determining  the  spa¬ 
tial  variability  of  the  wavefield  in  both  the  along-wave  and  the  cross-wave 
directions  was  the  primary  aim  of  the  data  analysis. 


PART  II:  DATA  COLLECTION 


Wave  Gages 

4.  Between  23  April  and  10  June  1983,  wave  data  were  collected  from  an 
array  of  five  Datawell  Waverider  wave  gages  moored  in  a  depth  of  17  m.  The 
array  was  positioned  about  2  km  offshore  of  the  FRF  (Figure  1).  The  dimen¬ 
sions  and  orientation  of  the  Waverider  array,  as  shown  in  Figure  2,  were  de¬ 
termined  after  placement  using  a  land-based  laser  surveying  instrument  during 
calm  wave  conditions.  All  the  buoys  in  the  array  were  tethered  to  mooring 
lines  of  the  same  scope  with  the  intent  that  any  displacement  due  to  a  unidi¬ 
rectional  current  would  result  in  a  translation  of  the  array  with  minimum  dis¬ 
tortion.  There  was  variation  in  the  array  dimensions  and  orientation  to  a 
limited  extent  due  to  wave-induced  displacements  of  the  buoys,  but  these  spa¬ 
tial  differences  averaged  out  over  the  period  of  data  collection. 

5.  Prior  to  deployment  all  the  Waveriders  were  calibrated  according  to 
the  manufacturer's  recommendations  using  a  Ferris  wheel  arrangement.  All  the 
gages  met  or  exceeded  the  calibration  specification  required  by  Datawell.  A 
posttest  check  of  the  calibration  curves  for  two  CERC  buoys  used  in  the  ex¬ 
periment  indicated  that  the  calibration  error  for  these  two  was  less  than 

4  percent. 


Data  Collected 

6.  The  data  set  was  comprised  of  38  data  collection  sessions  lasting 

80  min  each.  These  are  summarized  in  Table  1.  Data  were  obtained  for  periods 
of  wind-generated  seas,  swell,  and  for  cases  when  more  than  one  distinct  wave- 
train  was  present.  Each  session  involved  the  simultaneous  recording  of  wave 
data  from  all  gages  at  a  sampling  rate  of  four  points  per  second  resulting  in 
a  total  record  of  19,200  data  points  per  gage.  The  collection  sessions  (runs) 
were  substantially  longer  than  normal  because  it  was  realized  a  unique  data 
set  was  being  collected  and  that  the  longer  records  might  be  very  useful  for 
future  analyses.  The  raw  data  were  collected  at  the  FRF  on  a  Data  General 
Nova-4  computer  and  stored  on  magnetic  tape. 

7.  The  wave  direction  for  each  recording  session  was  routinely  obtained 
using  the  X-band  radar  at  the  FRF,  supplemented  by  visual  observations. 
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Figure  2.  Waverider  wave  gage  array 

Mattie  and  Harris  (1979)  report  the  direction  measuring  capability  of  the  *  cF 
radar  to  be  accurate  to  within  ±4  percent.  In  cases  with  two  distinct  wave- 
trains  approaching  from  different  directions,  the  directions  of  both  primary 
and  secondary  wavetrains  were  recorded.  Pertinent  weather  data  were  noted 
during  each  recording  session,  and  aerial  photographs  of  the  wave  field  were 
obtained  several  times  over  the  duration  of  the  experiment. 

8.  The  most  southerly  Waverider  buoy  (see  Figure  2)  produced  erratic 
results  from  the  onset  of  the  experiment  and  was  replaced  as  soon  as  possible 
Consequently,  the  data  from  the  first  22  collection  runs  do  not  include  data 
from  this  buoy.  Runs  after  run  22  include  results  from  the  most  southerly 
gage  in  the  array. 


PART  III:  DATA  ANALYSIS 


9.  Nine  of  the  data  collection  runs  were  selected  for  analysis,  as  in¬ 
dicated  in  Table  1,  and  represent  a  fairly  broad  range  of  conditions.  Further 
details  of  the  analysis  subset  are  provided  in  Table  2. 

10.  The  analysis  for  each  run  consisted  of  determining  the  energy  den¬ 
sity  spectrum  for  each  gage  and  performing  cross-spectral  analysis  to  deter¬ 
mine  wave  coherence  and  phase  functions  from  data  gathered  between  all  possi¬ 
ble  gage  pairings.  The  coherence  function  between  two  spatially  separated 
time  series  of  ocean  waves  measures  the  extent  to  which  the  time  series  at  one 
position  can  be  predicted  from  the  wave  field  at  the  other  position,  assuming 
a  linear  dependence  exists  between  the  two.  By  definition  the  coherence  func¬ 
tion  is  given  as 


Y _ (f)  = 


G  (f) 

xy 


G  (f)G  (f) 
xx  yy 


where 


G  (f)  =  cross-spectral  density  function  between  gages  X  and  Y 
xy 

G^^lf)  =  autospectral  density  function  of  gage  X 
G^(f)  =  autospectral  density  function  of  gage  Y 
The  coherence  values  are  bounded  between  zero  (no  relation  between  the  posi¬ 
tions)  and  a  coherence  of  one  (complete  linear  dependence).  The  coherence 
function  rapidly  decreases  as  the  distance  between  the  two  spatial  positions 
increases.  This  is  one  reason  that  ocean  waves  are  more  often  described  in 
terms  of  energy  spectra  with  parameters  relating  the  spectra  between  locations. 

11.  The  first  five  analysis  runs  resulted  in  six  possible  pairings  of 
gages  per  run  due  to  the  malfunctioning  of  the  most  southerly  Waverider.  The 
remaining  four  runs  each  contain  10  possible  pairings.  In  all  instances  only 
every  other  data  point  in  the  first  2048  sec  of  the  recorded  time  series  was 
used  in  the  data  analysis  subset.  Thus  the  analysis  subset  contains  4096 
points  sampled  at  0.5-sec  time  intervals.  This  provides  a  sufficient  number 
of  points  to  insure  accurate  representation  of  the  spectra  in  the  frequency 
range  0.05  to  1.0  Hz  (1-  to  20-sec  wave  periods)  without  becoming  overly  bur¬ 
densome  for  analysis.  There  is  evidence  that  coherence  estimates  are  not 
appreciably  improved  by  increasing  the  length  of  the  time  series  beyond  2048 
data  points  (Briggs  1981).  In  fact,  a  long  data  set  could  introduce  questions 
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about  stationarity  of  the  wave  field  during  the  recording  interval. 

12.  The  data  were  processed  with  a  Fast  Fourier  Transform  (FFT)  rou¬ 
tine.  The  resolution  in  the  frequency  domain  is  0.0004883  Hz.  Since  frequen¬ 
cies  higher  than  0.25  Hz  (4  sec)  are  considered  to  be  of  little  practical  im¬ 
portance  and  since  the  coherences  at  these  high  frequencies  are  expected  to  be 
very  small,  attention  was  concentrated  on  the  frequencies  ranging  from  0  to 
0.25  Hz.  The  analysis  produced  spectral  values  at  512  discrete  frequencies  in 
this  frequency  range. 

13.  Gages  were  paired  for  cross-spectral  analysis,  and  cross-spectral 
energy  densities  were  computed.  A  Tukey  filter  was  applied  to  the  cross¬ 
spectrum  Fourier  coefficients  to  increase  the  statistical  stability  of  the 
results  with  only  a  minor  loss  in  resolution.  Coherence  and  phase  were  com¬ 
puted  from  the  smoothed  cross-spectral  values. 

14.  It  would  have  been  possible  to  time  shift  one  time  series  with  re¬ 
spect  to  the  other  time  series  of  the  wave  gage  pairing  in  an  attempt  to  im¬ 
prove  the  coherence.  This  was  not  done  for  several  reasons.  First,  there  is 
a  difficulty  in  selecting  an  appropriate  time  shift  for  each  pairing.  A  pos¬ 
sible  candidate  would  be  the  time  it  takes  the  wave  associated  with  the  peak 
of  the  spectrum  to  traverse  the  distance  between  gages.  However,  this  lag 
time  will  vary  for  each  gage  and  doesn't  fully  compensate  for  all  the  waves 
traveling  at  different  speeds.  Second,  the  amount  of  wave  record  not  over¬ 
lapped  in  time  (based  on  the  travel  time  of  the  peak  period  wave)  never  ex¬ 
ceeded  1.5  percent  of  the  total  record  and  most  often  this  percentage  was  much 
less.  While  not  actually  tested,  it  was  believed  that  time  shifting  one  rec¬ 
ord  with  respect  to  another  would  not  significantly  alter  the  coherence,  and 
it  would  have  the  disadvantage  of  introducing  an  additional  questionable  as¬ 
sumption  into  the  analysis. 

15.  Coherence  values  greater  than  0.7  (coherence-squared  values  greater 
than  about  0.5)  are  considered  to  indicate  a  meaningful  coherence  between  sig¬ 
nals,  although  admittedly  this  is  a  rather  arbitrary  demarcation.  The  anal¬ 
yses  were  searched  for  coherence  values  greater  than  0.7.  Fifty  percent  of 
the  cross-spectral  analysis  indicated  such  coherences  in  at  least  part  of  the 
frequency  range.  The  cases  which  did  not  have  any  coherence  values  above  0.7 
represented  pairing  of  the  gages  which  had  large  spatial  separation  (partic¬ 
ularly  in  the  along-crest  direction)  for  the  incident  wave  direction.  For 
all  cases  with  reasonably  high  coherence,  the  coherence  function  and  the 


autospectral  energy  density  functions  were  plotted,  as  illustrated  in  Fig¬ 
ure  3.  The  entire  set  of  plots  is  given  in  Appendix  A.  For  plotting  purposes 
the  discrete  frequency  intervals  were  combir**d  into  100  bands  by  averaging 
every  five  values.  This  had  only  a  minor  smoothing  effect  on  the  plotted 
results . 
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Figure  3.  Example  of  the  Appendix  A  wave  spectra  and  coherence  plots 
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PART  IV:  RESULTS 

Classification  in  Terms  of  Incident  Wave  Direction 

16.  For  every  incident  wave  direction,  the  vector  distance  between  each 
wave  gage  pairing  can  be  resolved  into  a  component  parallel  to  the  wave  crests 
and  a  component  perpendicular  to  the  crests  as  illustrated  by  Figure  4.  This 
breakdown  provides  a  means  of  assessing  spatial  variability  in  the  wave  field. 
The  "along-crest"  distance  has  been  labeled  as  Y,  and  the  "down-crest"  dis¬ 
tance  as  X  on  the  plots  of  coherence  in  Appendix  A.  Both  X  and  Y  are  ex¬ 
pressed  in  meters.  Values  of  X  and  Y  are  given  for  both  the  primary  and  the 
secondary  directions  when  two  distinct  wave  trains  were  present  (see  Fig¬ 
ure  4).  HMO  on  the  plots  is  equal  to  four  times  the  standard  deviation  of  sea 
surface  elevations  in  meters,  but  is  commonly  referred  to  as  significant  wave 
height.  The  gage  numbers  of  the  pairing  used  for  the  plot  refer  to  the  Wave- 
rider  numbering  as  shown  in  Figure  2. 

Coherence  Trends 

17.  Several  qualitative  trends  are  observed  in  the  plots  in  Appendix  A. 
The  coherence  function  is  weakly,  if  at  all,  dependent  upon  the  shape  of  the 
spectrum  or  the  amount  of  energy  contained  in  the  spectrum.  For  example,  the 
two  plots  given  in  Figure  5  show  coherence  functions  which  are  quite  similar 
even  though  the  spectra  between  run  3  and  run  7  are  considerably  different. 

Run  3  represents  a  narrow  spectrum  while  run  7  is  a  broad  spectrum  with  the 
energy  spread  over  a  wide  range  of  frequencies. 

18.  The  coherence  function  does  not  appear  to  be  influenced  greatly  by 
the  width  of  the  spectral  peak  or  by  the  frequency  at  which  it  is  located. 

The  region  of  good  coherence  for  run  3  in  Figure  5  is  considerably  wider  than 
the  peak  of  the  spectrum.  Run  6  (Figure  6)  represents  spectra  with  peak  fre¬ 
quencies  higher  than  those  in  Figure  5  but  with  a  very  similar  coherence 


function. 


Spatial  Variability 


19.  The  most  apparent  influence  on  the  coherence  function  is  under¬ 
standably  the  spatial  separation  of  the  buoys.  The  coherence  plots  in 
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Appendix  A  have  been  arranged  in  order  of  increasing  along-crest  (Y)  separa¬ 
tion  distances.  The  trend  is  a  decrease  in  coherence,  particularly  in  the 
higher  frequencies,  as  the  along-crest  separation  distance  increases,  some¬ 
times  regardless  of  the  down-crest  (X)  separation.  While  the  down-crest  dis¬ 
tance  has  some  influence,  the  overall  variation  in  coherence  function  is  more 
closely  tied  to  the  along-crest  distance  than  to  any  other  single  parameter. 

A  convenient  way  to  represent  this  "short-crestedness"  is  to  nondimensionalize 
the  along-crest  (Y)  and  the  down-crest  (X)  distances  by  the  wavelength  (L) 
associated  with  each  frequency  in  17-m  water  depth.  All  discrete  values  of 
coherence  can  then  be  plotted  against  the  separation  distance  in  terms  of  rel¬ 
ative  wavelength. 

20.  To  illustrate  the  importance  of  the  along-crest  separation  versus 
the  down-crest  separation,  plots  were  made  using  all  discrete  coherence  values 
between  0.078  and  0.16  Hz,  which  covers  the  range  of  meaningful  coherence  and 
the  principal  energy-containing  frequencies  in  nearly  all  cases.  Figure  7 
gives  coherence  versus  the  nondimens ional  along-crest  distance  for  all  values 
except  the  multiwavetrain  cases,  and  Figure  8  contains  the  same  coherences  as 
a  function  of  nondimensional  down-crest  distance.  A  tendency  for  coherence  to 
decrease  more  rapidly  with  increasing  along-crest  separation  than  with  in¬ 
creasing  down-crest  separation  is  evident  in  the  figures.  The  line  shown  in 
Figure  7  could  represent  an  approximate  upper  limit  boundary  to  the  coherence 
data.  Up  to  values  of  Y/L  around  1.5,  the  curve  is  probably  a  good  estimate 
of  the  best  possible  coherence  which  could  be  obtained  between  two  gages 
aligned  parallel  to  the  wave  crests.  A  similar  upper  limit  curve  is  drawn  in 
Figure  8,  but  the  plotted  values  of  coherence  are  more  evenly  distributed  be¬ 
tween  zero  and  one  than  those  in  Figure  7.  Thus  the  upper  limit  line  in  Fig¬ 
ure  8  is  less  well  defined.  These  upper  limit  curves  should  be  considered 
only  in  conjunction  with  spatial  separations  on  the  order  of  30  m  or  more. 

21.  Attempts  were  made  to  examine  the  data  as  a  function  of  the  magni¬ 
tude  of  the  separation  vector  and  the  angle  of  orientation  of  the  vector  with 
respect  to  the  incident  waves.  Qualitatively,  the  data  followed  the  expected 
trends  but  the  scatter  was  far  too  wide  to  justify  any  quantitative  treatment. 

Coherence  Contours 

22.  A  nondimensional  plot  of  coherence  contours  is  given  in  Figure  9. 
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Figure  7.  Coherence  versus  nondimensional  along-crest  separation 


down-crest  separation 


The  horizontal  axis  gives  the  along-crest  separation  distance  (Y)  divided  by 
the  wavelength  (L)  associated  with  each  frequency  (linear  theory),  while  the 
vertical  axis  gives  the  nondimens ional  down-crest  separation  (X/L) .  The  nu¬ 
merical  values  on  the  plot  are  the  average  of  all  discrete  values  of  coherence 
which  fell  into  0.2  by  0.2  squares  centered  on  the  decimal  points  on  the  plot. 
The  underlined  values  represent  the  squares  containing  more  than  100  points  in 
the  average,  and  the  dashed  sections  of  the  contours  are  surmised.  The  number 
of  coherence  values  used  in  obtaining  each  average  was  subjectively  considered 
when  sketching  the  contours.  The  multiwavetrain  data  were  not  included  in 
this  formulation;  thus,  the  contours  are  intended  for  a  wave  field  approaching 
from  a  single  central  direction. 

23.  The  contours  are  a  crude  estimate  at  best  since  they  consider  only 
the  spatial  separation  parameters  of  a  complex  process.  However,  they  do  pro¬ 
vide  a  means  of  getting  a  first  estimate  of  the  coherence  function  between  two 
gage  locations.  The  Y-  and  X-axis  intercept  values  of  the  contour  lines  in 
Figure  9  are  somewhat  less  than  values  obtained  from  the  upper  limit  lines  in 
Figures  7  and  8  because  the  intercept  values  are  meant  to  represent  average 
values  instead  of  upper  limit  values. 

24.  To  illustrate  the  use  of  the  contours,  four  coherence  plots  from 
Appendix  A  were  chosen,  and  an  estimate  of  coherence  was  made  using  the  con¬ 
tours  of  Figure  9.  These  results  are  shown  as  dashed  lines  in  Figure  10.  The 
estimates  are  fair  in  the  region  of  highest  coherence,  but  fail  to  give  satis¬ 
factory  reproduction  at  the  higher  frequencies.  It  should  again  be  noted  that 
only  the  frequency  range  0.078  to  0.16  Hz  was  considered  in  preparing  the 
contours . 

25.  It  must  be  emphasized  that  the  results  which  were  contoured  in  Fig¬ 
ure  9  were  obtained  in  a  water  depth  of  17  m.  The  contours  may  be  useful  at 
other  shallow  water  depths  because  of  the  nondimensionalization  by  wavelength; 
however,  more  data  are  needed  to  support  this  generalization. 

Analytical  Phase  Shift 

26.  An  analytical  experiment  aimed  at  obtaining  improved  coherence  be¬ 
tween  two  wave  gages  was  performed  on  three  gage  pairs.  These  pairings  were 
selected  because  they  came  the  closest  to  having  one  gage  situated  directly 
down-crest  of  the  other  gage  in  the  pair.  In  other  words,  the  "along-crest" 
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distance  was  minimal.  The  experiment  consisted  of  deterministically  trans¬ 
lating  the  upwind  wave  gage  results  to  the  downwind  gage  position  using  the 
linear  dispersion  relation.  It  was  hoped  this  would  compensate  for  the  fact 
that  waves  at  different  frequencies  travel  at  different  speeds.  A  coherence 
function  could  then  be  determined  between  the  downwind  gage  and  the  phase- 
shifted  results  from  the  upwind  gage.  Then  a  comparison  could  be  made  with 
the  original  coherence  results  to  see  if  any  improvement  in  coherence  could 
be  realized.  A  successful  experiment  would  allow  an  improved  estimate  of  the 
wavetrain  at  a  location  based  on  wave  data  collected  a  short  distance  away. 

27.  The  phase  shifting  was  accomplished  by  modifying  the  upwind  gage 
Fourier  coefficients  obtained  from  the  FFT.  For  each  discrete  frequency  w 
a  pair  of  Fourier  coefficients,  Pj(iu)  and  Qj(u>)  ,  are  obtained.  They  re¬ 
late  to  the  wave  spectrum  through  the  following  equations: 

i0  (w) 

Fj  (in)  =  P1  (in)  +  iQ  j  (in)  =  jF^m)!  e  (2) 

where 

l  2  2 

Fj(in)  =  .JP^On)  +  Q^in)  =  discrete  component  of  the  (3) 

amplitude  spectrum 


QjOju) 


tan 


-1 


Qj  (u») 

Pj(U)) 


=  discrete  component  of  the 
phase  spectrum 


(4) 


It  was  assumed  for  the  experiment  that  the  phase  angle  at  the  downwind  posi¬ 
tion  would  be 


02(w)  =  (in)  +  A0(u>) 
Since  the  phase  angle  can  be  expressed  as 

0(u>)  =  int  +  kx 


(5) 


(6) 


with 


in  =  circular  frequency 
t  =  time 

x  =  downwind  spatial  distance 


k  =  wave  number  associated  with  the  frequency  u>  using  linear  theory 
the  spatial  correction  to  the  phase  angle  for  t  =  constant  becomes 


A6(u>)  _ 
Ax 


or 


A8(u>)  =  kAx 


(7) 


where  Ax  is  the  spatial  separation  of  the  wave  gages. 

28.  To  implement  this  phase  it  is  necessary  at  each  frequency  to  deter¬ 
mine  a  new  set  of  Fourier  coefficients  giving  the  new  phase  angle  without 
changing  the  amplitude,  i.e., 

02(u>)  =  0j(u>)  +  kAx  =  tan 

and 


=  iF^uOl  (9) 

Thus  the  new  coefficients  become 


P2(u>)  = 

(F1M|  cos  | 

fej  <■••!)  +  kAxJ 

j  OO) 

Q2(w)  = 

j  Fj  (u»)  |  sin  | 

^(u)  +  kAxJ 

|  (ID 

29.  Figures  11,  12,  and  13  display  the  results  of  this  experiment  with 
the  analytically  shifted  coherence  given  by  the  dashed  line.  Obviously,  the 
concept  is  not  a  workable  scheme  in  reality.  A  possible  source  of  the  diffi¬ 
culty  is  the  use  of  the  linear  dispersion  relationship. 

30.  A  more  probable  cause  for  the  failure  of  the  above-mentioned  phase 
modification  lies  with  the  phase  spectra  resulting  from  the  FFT .  Although  the 
phase  spectra  produced  by  the  FFT  can  be  combined  with  the  corresponding  am¬ 
plitude  spectra  to  recover  the  original  time  series,  it  has  been  demonstrated 
by  Thompson  (1982)  that  phase  spectra  from  FFT  analyses  may  not  be  a  true 
representation  of  the  phases  associated  with  the  physical  process  being  ana¬ 
lyzed.  Phase  spectra  from  FFT  analyses  are  characterized  by  large  variations 
between  adjacent  frequencies.  Thompson's  research  showed  that  these  fluctua¬ 
tions  can  be  induced  by  limitations  of  the  FFT  analysis  and  are  thus  partial 
artifacts  of  the  procedure  itself.  The  fluctuations  appear  to  be  a  conse¬ 
quence  of  representing  the  physical  process  by  a  finite  number  of  discrete 
frequency  components.  Therefore,  modification  of  an  FFT-produced  phase 
spectra  is  likely  to  be  unsuccessful. 
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PART  V:  SUMMARY  AND  RECOMMENDATIONS 


31.  The  primary  objective  of  this  experiment  was  to  measure  the  spatial 
variability  in  the  wave  field  under  a  variety  of  wave  conditions  so  that  find¬ 
ings  could  assist  future  ship  mooring  testing  by  the  NCEL.  This  report  has 
detailed  the  experiment  and  the  results  obtained  directly  from  the  measure¬ 
ments  and  those  obtained  through  manipulation  of  the  data  set.  Reference  is 
made  to  the  entire  set  of  plots  in  Appendix  A.  The  results  are  presented  in 
terms  of  the  wave  coherence  function. 

32.  Qualitatively,  the  coherence  function  appears  relatively  indepen¬ 
dent  of  the  shape  of  the  spectrum,  the  spectral  energy  content,  and  the  loca¬ 
tion  of  the  spectral  peak  frequency  for  the  range  of  conditions  examined.  The 
finite  water  depth  at  the  site  appears  to  render  the  lower  frequency  waves 
less  dispersive  and,  hence,  more  coherent  so  long  as  sufficient  energy  is 
present  at  those  frequencies. 

33.  The  most  apparent  influence  on  the  coherence  is  the  spatial  separa¬ 
tion,  particularly  in  the  along-crest  direction.  While  too  much  scatter  is 
present  to  properly  quantify  the  effect  of  separation,  an  upper  limit  is 
given.  It  is  also  possible  to  construct  a  set  of  nondimensionalized  coherence 
contours  for  the  purpose  of  estimating  the  coherence  function  in  shallow 
water.  The  contours  provide  a  fair  estimate  in  the  region  of  good  coherence 
when  compared  with  actual  records,  but  they  overestimate  the  coherence  as  the 
frequency  increases  beyond  about  0.13  Hz. 

34.  The  experiment  provided  some  insight  into  the  variability  of  the 
nearshore  wave  field,  and  from  this  knowledge  it  is  possible  to  make  some  op¬ 
erational  recommendations  to  the  NCEL  with  regard  to  their  future  experiments. 

35.  The  contours  developed  in  Part  IV  will  be  helpful  in  planning  ship 
mooring  experiments,  but  it  is  evident  that  the  NCEL  should  also  plan  to  col¬ 
lect  onsite  data  to  determine  wave  coherence.  The  results  indicate  that  wave 
coherence  rapidly  decreases  as  spatial  separation  increases,  particularly  in 
the  along-crest  direction. 

36.  As  suggested  by  the  NCEL  project  engineer,  a  good  experimental  pro¬ 
cedure  would  be  to  place  one  Waverider  buoy  a  distance  d^  from  the  ship  and 
a  second  buoy  a  distance  of  d^  from  the  first  buoy.  These  buoys  should  be 
in  line  with  the  ship  and  oriented  directly  "up-crest"  from  the  vessel.  In 
this  manner  wave  coherence  will  be  maximized  since  the  along-crest  separation 


is  minimal,  and  the  coherence  function  found  between  the  two  gages  can  be  as¬ 
sumed  to  be  reasonably  close  to  the  coherence  between  the  waves  at  the  buoy 
nearest  to  the  ship  and  those  causing  the  ship's  motion. 

37.  A  certain  amount  of  variation  in  incoming  wave  direction  can  be 
tolerated  without  significantly  reducing  wave  coherence,  but  it  will  probably 
be  necessary  to  realign  the  Waveriders  several  times  over  the  course  of  an  ex¬ 
periment  as  wave  direction  changes.  The  contours  in  Figure  9  can  be  used  to 
estimate  allowable  limits  of  wave  approach  angle  variation  before  the  buoys 
should  be  moved.  To  aid  in  the  Waverider  relocation,  it  may  be  advisable  to 
have  an  array  of  vacant  buoy  mooring  installed  at  the  onset  of  the  experiment. 
Since  the  mooring  array  spacings  and  orientation  would  have  already  been  de¬ 
termined,  the  relocation  of  buoys  could  proceed  more  smoothly,  particularly  in 
rougher  weather.  The  mooring  array  configuration  can  be  optimized  by  the  de¬ 
termination  of  "typical"  seasonal  wave  directions  for  the  particular  site. 

38.  It  is  highly  advisable  that  the  NCEL  have  at  least  one  backup  Wave- 
rider  gage  moored  at  the  site  to  insure  a  relative  continuity  in  the  data  col¬ 
lection.  The  ideal  situation  would  be  to  have  two  backup  gages  moored  in  the 
array  in  such  a  fashion  as  to  become  the  primary  gages  when  the  waves  approach 
perpendicularly  to  the  gages'  alignment. 

39.  There  is  no  clear  difference  in  wave  coherence  between  conditions 
representing  swell  and  storm  seas  in  the  data  analyzed.  Swell  waves  can  be 
expected  to  be  longer  crested  and  less  dispersive  and  thus  less  sensitive  to 
spatial  separation,  but  this  was  not  clearly  evident  in  the  data.  This  obser¬ 
vation  may  be  a  consequence  of  the  finite  depth  influence,  and  the  difference 
between  sea  and  swell  might  become  more  apparent  in  deeper  water. 

40.  All  of  the  results  were  obtained  in  a  water  depth  of  17  m.  There¬ 
fore,  it  is  reasonable  to  expect  some  signature  of  the  depth  in  the  coherence 
function.  The  high  coherence  values  observed  in  the  frequency  range  0.08  to 
0.1  Hz  are  quite  possibly  related  to  the  fact  that  waves  at  these  frequencies 
have  been  rendered  less  dispersive  by  the  finite  depth.  Although  it  is  impos¬ 
sible  to  prove  with  only  the  existing  data  set,  it  is  reasonable  to  hypothe¬ 
size  that  the  high  coherences  at  relatively  low  frequencies  are  a  result  of 
the  waves  progressing  with  little  alteration  in  form.  According  to  linear 
wave  theory,  the  ratio  of  phase  velocity  to  group  velocity  for  a  0.1-Hz 
(10-sec)  wave  in  17-m  depth  is  1.26.  The  observed  decrease  in  coherence 
values  at  frequencies  less  than  about  0.07  to  0.08  Hz  may  well  be  due  to  the 


minuscule  amounts  of  energy  at  these  low  frequencies. 

41.  The  data  analyzed  in  this  study  provide  strong  evidence  of  higher 
coherences  at  frequencies  less  than  about  0.125  Hz  (8-sec  period)  than  at 
higher  frequencies.  This  trend  appears  unrelated  to  the  form  of  the  autospec¬ 
trum.  Since  one  intention  of  mooring  experiments  is  to  have  high  coherence  at 
energy-containing  frequencies  affecting  ship  motion,  wave  fields  with  energy 
primarily  at  frequencies  lower  than  0.125  Hz  should  be  preferred  for  the 
experiment. 

42.  The  mooring  lines  for  the  Waveriders  should  have  the  same  scope  to 
minimize  the  variation  in  buoy  separation  distance.  The  wave-induced  buoy 
displacements  are  not  thought  to  be  a  factor  in  the  large  fluctuations  in  the 
coherence  function.  This  is  supported  by  Kuo,  Mitsuyasu,  and  Masuda  (1979a, 
1979b)  in  which  the  same  coherence  fluctuations  appear  between  wave  gages 
which  have  been  rigidly  fixed  in  position.  Change  in  separation  distance  be¬ 
tween  the  moored  ship  and  the  nearest  Waverider  will  depend  on  the  freedom  of 
movement  in  the  ship's  mooring  system  and  the  ship's  dynamic  response  to  the 
mooring  line  forces. 

43.  Sea  states  having  multiple  wavetrains  approaching  from  different 
directions  exhibit  wave  coherence  functions  which  deteriorate  more  rapidly 
than  sea  states  with  single  wavetrains  and  a  single  central  direction.  When 
recording  multiple  wavetrains,  it  is  best  to  align  with  the  direction  contain¬ 
ing  the  most  energy.  However,  the  presence  of  multiple  wavetrains  will  cer¬ 
tainly  complicate  data  analysis  and  interpretation. 

44.  The  following  recommendations  are  offered  with  regard  to  future 
NCEL  ship  mooring  experiments . 

a.  Align  the  two  Waverider  buoys  parallel  with  the  central  inci¬ 
dent  wave  direction. 

b.  Place  additional  Waverider  moorings  in  the  array  to  facilitate 
buoy  relocation. 

c.  Relocate  Waveriders  when  feasible  to  maintain  the  desired 
alignment  with  the  wave  field. 

d.  Maintain  equal  scope  on  all  the  wave  gage  moorings. 

e.  Moor  at  least  one  (more  if  possible)  backup  Waverider  at  the 
experiment  site. 


f.  If  a  choice  exists,  avoid  collecting  experimental  data  during 
wave  conditions  with  energy  concentrated  at  the  higher  frequen¬ 
cies  (greater  than  0.13  Hz)  and  during  sea  states  with  multiple 
wavetrains  approaching  from  different  directions. 
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APPENDIX  A:  COHERENCE  PLOTS 


A2 


ENERGY  DENSITY  VS  FREQUENCY 


RUN  •  1 
1ST  GAGE  •  1 
2ND  GAGE  •  4 


O.M.OI  0.02  0.03  0.04  0.05  0.06  0.07  0.00  0.09  0.10  0.11  0.12  0.13  0.14  0.15  0.  IS  0.17  0.10  0.19  0.20  0.21  0.22  0.2J  0.24  0, 

FREQUENCY- -HZ 


COHERENCE  VS  FREQUENCY 


huo  -  0.5 
X  -  67.4 
Y  -  151.2 
XSEC  *  163.4 
YSEC  -  25.9 


O.OS.OI  0.02  0.03  0.04  0.05  0.00  0.07  0.00  0.00  0.10  0.1 1  0.12  0.13  0.14  0.15  0.10  0.17  0.10  0.19  0.20  0.21  0.22  0.23  0.24  0 

FREQUENCY— HZ _ 


Figure  A2.  Analysis  run  1  of  gages  1  and  4 


V'A.  '  W  V  V  ’.•  ' 


-■.jfVfi'O  r 


O.M.OI  0.03  0.01  0.01  0.05  0.06  0.07  0.06  0.00  0.10  0.11  0.12  0.13  0.11  0.15  0.16  0.17  0.16  0.16  0.30  0.21  0.23  0.33  0.21  0, 

_ frequency — HZ  _ 


Figure  A4.  Analysis  run  2  of  gages  0  and  4 


ENERGY  DENSITY  VS  FREQUENCY 


RUN  •  2 
1ST  GflGC  •  0 
2ND  GflGE  •  l 


o.oa.oi  0.02  O.OS  0.04  O.OS  0.00  0.07  0.00  0.00  0.10  0.1 1  0. 12  0.  IS  0.14  0.  IS  0.10  0.17  0.10  O.IS  0.20  0.21  0.22  0.23  0.24  0 


run  2  of  gages  0  and  1 


RUN  «  2 
1ST  GflGE  •  1 
2ND  GAGE  •  4 


ENERGY  DENSITY  VS  FREQUENCY 


O.OB.OI  0.02  0.03  0.09  0.05  0.06  0.07  0.08  0.09  0.10  0.11  0.12  0.13  0.19  O.IS  0.16  0.17  0.18  0. 19  0.20  0.21  0.22  0.23  0.29  0, 

FREQUENCY — HZ 


HHO  -  0.8 
X  -  80.2 
Y  •  144. 8 


COHERENCE  VS  FREQUENCY 


•  v 


O.OB.OI  0.02  0.03  0.09  0.05  0.08  0.07  0.08  0.09  0.10  O.lt  0.12  0.13  0.19  0.15  0.16  0.17  0.18  0.19  0.20  0.21  0.22  0.23  0.29  0. 

_ _ FREQUENCY— HZ _ 


Figure  A6.  Analysis  run  2  of  gages  1  and  4 


Figure  A10.  Analysis  run  3  of  gages  0  and  1 


All 


RUN  a  3 
1ST  GflGE  a  3 
2ND  GflGE  a  4 


ENERGY  DENSITY  VS  FREQUENCY 


A12 
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ENERGY  DENSITY  VS  FREQUENCY 


RUN  •  4 
1ST  Gfl6E  •  1 
2ND  GflGE  «  3 


0.00.01  0.02  0.01  0.04  0.05  0.00  0.0?  0.00  0.09  0.10  0.11  0.12  0.13  0.14  0.15  0.10  0.17  0.19  0.11  0.20  0.21  0.22  0.23  0.24  0 

fr 


run  4  of  gages  1  and  3 


D_r. 


Figure  A13.  Analysis  run  4  of  gages  0  and  4 


ENERGY  DENSITY  VS  FREQUENCY 


RUN  •  4 
1ST  GflGE  •  0 
2ND  GflGE  »  1 


o.aa.01  o.o2  o.as  o.m  0.05  0.06  0.07  0.0a  0.09  0.10  o.tt  0.12  0.13  0.14  0.1s  0.16  0.17  0. is  0.19  0.20  0.21  0.22  0.23  0.24  0 

FREQUENCY— HZ 


COHERENCE  VS  FREQUENCY 


HMO  -  1.6 
X  -  55.5 
Y  -  72.8 


O.M.01  0.02  0.09  0.01  0.05  0.00  0.07  0.00  0.09  0.10  O.II  0.12  0.13  0.14  0.19  0.(6  0.17  0.10  0.11  0.20  0.21  0.22  0.23  0.24  0. 

_ _ FREQUENCY— HZ _ 

Figure  A14.  Analysis  run  4  of  gages  0  and  1 


art  km  oo*i  wo  wo  tro 


0.25  0.50  0.T5  1.00 


ENERGY  DENSITY  VS  FREQUENCY 


RUN  •  5 
1ST  GAGE  •  0 
2ND  GAGE  •  3 


0.00. 01  0.02  0.03  0.04  0.05  0.06  0.07  0.08  0.09  0.10  O.lt  0.12  0.13  0.14  0.15  0.16  0. 17  0.16  0.19  0.20  0.21  0.22  0.23  0.24  0 

FREQUENCY— HZ 


COHERENCE  VS  FREQUENCY 


HO  -  |.I 

-  2.7 
•  137.2 
SEC  -  135.6 
SEC  -  21.0 


0.oa.0i  0.02  0.03  0.01  0.05  0.06  0.07  0.06  0.09  0.10  O.lt  0.12  0.13  0.14  0.15  0.16  0.17  0.16  0.19  0.20  0.21  0.22  0.23  0.24  0 

_ TRE 


Figure  A17.  Analysis  run  5  of  gages  0  and  3 


8 


O.M.OI  0.03  0.03  0.04  0.05  0.00  0.07  0.00  0.01  0. 10  0. 1 1  0. 13  0. 13  0. 14  0.  IS  0. 10  0. 17  0. 10  0.  IS  0.30  0.31  0.33  0.33  0.24  0 


run  6  of  gages  0  and  4 


RUN  •  6 
1ST  GflGE  •  0 
2ND  GflGE  •  1 


ENERGY  OENSITY  VS  FREQUENCY 


COHERENCE  VS  FREQUENCY 


0.25  0.50  0.75  1.00 


ENERGY  DENSITY  VS  FREQUENCY 


COHERENCE  VS  FREQUENCY 


A24 


0.15  0.50  0.79  1.00 


ENERGY  DENSITY  VS  FREQUENCY 


0.M.01  0.02  0.01  O.M  0.05  0.06  0.07  0.00  0.09  0.10  O.lt  0.12  0.13  0.11  0.15  0.16  0.17  0.10  0.19  0.20  0.21  0.22  0.23  0.21  0. 


FREQUENCY- -HZ 


COHERENCE  VS  FREQUENCY 


Figure  A24.  Analysis  run  7  of  gages  1  and  4 


0.00  0.M  O.SO  O.n  1.00 


ENERGY  OENSITY  VS  FREQUENCY 


COHERENCE  VS  FREQUENCY 


S8 


ENERGY  DENSITY  VS  FREQUENCY 


RUN  •  8 
1ST  6RGE  •  2 
2ND  GflGE  •  A 


0.00.01  0.02  o.aa  0.04  o.os  o.t*  a.07  o.oo  0.09  o.io  o.tt  0.12  0.13  0.14  0.1s  o.io  0.17  o.io  o.ig  0.20  0.21  0.22  0.2s  0.24  0 

FREQUENCY— HZ 


COHERENCE  VS  FREQUENCY 


HW  -  1.1 
X  -  17.7 

Y  -  25.0 
XSEC  -  2.7 
YSEC  -  30.5 


01  0.02  0.09  0.04  0.05  0.00  0.07  0.00  0.00  0. 10  0. It  0. 12  0. 13  0. 14  0. 15  0. IS  0. 17  0. 10  0. 10  0.20  0.21  0.22  0.23  0.24  0 


Figure  A26.  Analysis  run  8  of  gages  2  and  4 


ENERGY  DENSITY  VS  FREQUENCY 


RUN  •  8 
1ST  6R6E  •  0 
2ND  GflGE  •  I 


O.OS.OI  0.02  0.03  0.M  0.05  0.00  0.07  0.00  0.00  0.10  O.il  0.13  0.13  0.14  0.15  0.10  0.17  0.10  0.15  0.20  0.21  0.22  0.23  0.24  0 

_ _ FREQUENCY— HZ 


Figure  A27.  Analysis  run  8  of  gages  0  and  1 


ENERGY  DENSITY  VS  FREQUENCY 


RUN  •  8 
1ST  GflGE  •  0 
2ND  GflGE  •  4 


0.09.01  0.02  0.03  0.M  0.05  O.OB  0.07  0.00  0.09  0. 10  0. 1 1  0. 12  0. 13  0. 14  0.  IS  0.  >6  0. 17  0. 19  0. 19  0.20  0.2!  0.22  0.23  0.24  OJ 

FREQUENCY--HZ 


COHERENCE  VS  FREQUENCY 


HNO  -  1.1 
X  -  42.7 
Y  -  63.1 
XSEC  *  68.6 
YSEC  -  33.5 


9.09.01  0.02  0.09  0.04  0.05  0.00  0.07  0.00  0.00  0. 10  0. 1 1  0. 12  0. 13  0. 14  0. 15  0. 10  0. 17  0. 10  0. 19  0.20  0.21  0.22  0.23  0.24  0, 

_ £SEQU6NCY-’HZ - 

Figure  A28.  Analysis  run  8  of  gages  0  and  4 
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A29 


RUN  «  9 
1ST  GAGE  »  0 
2ND  GAGE  >  1 


ENERGY  DENSITY  VS  FREQUENCY 


0.08.01  0.03  0.03  0.04  0.00  0.00  0.07  0.00  0.00  0.10  0.11  0.ia  0.11  0.14  0.10  0.10  0.17  0. II  0.10  0.30  0.31  0.33  0.33  0.34 


Figure  A29.  Analysis  run  9  of  gages  0  and  1 


Figure  A31.  Analysis  run  9  of  gages  0  and  4 


3.001.01  0.03  0.03  0.04  O.OS  0.36  0.07  0.06  0.09  0.10  0.11  0.13  0.13  0.14  0.15  0.16  0.17  0.16  0.16  0.30  0.31  0.33  0.33  0.34  0. 


0(11.25  0.50  0.75  1.00  1.25  1.50  1.75  2.00  2.25  2.50  2.75 


ENERGY  DENSITY  VS  FREQUENCY 


RUN  *  9 
1ST  GAGE  «  1 
2  NO  GAGE  •  2 


FREQUENCY-HZ 


COHERENCE  VS  FREQUENCY 


o.oa.oi  o.(U  0.03  o.w  o.os  o.os  o.o7  o.oe  o.oo  o.io  o.u  0.12  0.13  0.14  o.ts  a.  te  0.17  0.11  o.u  0.200.21  0.22  o.os  0.24  0 
_  _  PR 


s  run  9  of  gages  1  and  2 


0.02  0.03  0.01  3.05  0.06  0.07  0.06  0.09  0.  tO  0.11  0.12  0.13  0.11  O.lS  0.16  0.17  0.16  0.19  0.20  0.21  0.22  0.23  0.21  I 

FREQUENCY— HZ 


COHERENCE  VS  FREQUENCY 


Figure  A35.  Analysis  run  9  of  gages  1  and  3 
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